Bistability of vortex core dynamics in a single perpendicularly magnetized nano-disk 



ON 

o 
o 

(N 

>>: 



in 



C/5 



G. de Loubens^'B A. Riegler, 2 B. Pigeau, 1 F. Lochncr, 2 F. Boust, 3 K. Y. Guslienko, 4 H. Hurdcquint, 5 
L. W. Molenkamp, 2 G. Schmidt, 2 <Q A. N. Slavin, 6 V. S. Tiberkevich, 6 N. Vukadinovic, 7 and 0. Klein 1 

1 Service de Physique de I'Etat Condense (CNRS URA 2464), CEA Saclay, 91191 Gif-sur-Yvette, France 
2 Physikalisches Institut (EP3), Universitdt Wiirzburg, 97074 Wiirzburg, Germany 
ONERA, Chemin de la Huniere, 91761 Palaiseau, France 
^ Department of Materials Physics, The University of the Basque Country, 20080 San Sebastian, Spain 
Laboratoire de Physique des Solides, Universite Paris-Sud, 91405 Orsay, France 
''Department of Physics, Oakland University, Michigan 48309, USA 
7 Dassault Aviation, 78 quai Marcel Dassault, 92552 Saint-Cloud, France 
(Dated: May 5, 2009) 

Microwave spectroscopy of individual vortex-state magnetic nano-disks in a perpendicular bias 
magnetic field, H, is performed using a magnetic resonance force microscope (MRFM). It reveals 
the splitting induced by H on the gyrotropic frequency of the vortex core rotation related to the 
existence of the two stable polarities of the core. This splitting enables spectroscopic detection of 
the core polarity. The bistability extends up to a large negative (antiparallel to the core) value 
of the bias magnetic field H r , at which the core polarity is reversed. The difference between the 
frequencies of the two stable rotational modes corresponding to each core polarity is proportional to 
H and to the ratio of the disk thickness to its radius. Simple analytic theory in combination with 
micromagnetic simulations give quantitative description of the observed bistable dynamics. 
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Magnetic vortices are singular topological states found 
in the equilibrium magnetic configuration of sub-micron 
size ferromagnetic dots [J . In a certain range of dot 
aspect ratios (ratio j3 = t/R of the dot thickness t to its 
radius R) the equilibrium ground state of the static mag- 
netization consists of the curling in-plane magnetization 
and a nanometer size core of the out-of-plane magnetiza- 
tion at the dot center. The magnetization of the vortex 
core can point either up or down, both polarities p = ±1 
being degenerate at zero field. This bi-stable property of 
magnetic vortices, as well as the switching from one po- 
larity to the other, have been intensively studied in the 
past few years because of their possible applications in 
magnetic storage devices 0, 0, IE El • It has already been 
established : (i) that the lowest excitation mode of the 
vortex state is the gyrotropic mode corresponding to a 
rotation of the vortex core about the dot center, (ii) that 
the frequency of this mode is linearly proportional the 
dot aspect ratio (3 0] , and (iii) that the sense of gyration 
of the vortex core is determined by a right-hand rule to 
the core polarity [3] ■ 

In this Letter, we report that by using the exquisitely 
sensitive method of magnetic resonance force microscopy 
(MRFM) Q, we were able to observe bistability of the 
vortex core dynamics in a single magnetic disk subjected 
to a perpendicular bias magnetic field, that was varied in 
a wide range from positive (parallel to the vortex core) 
to negative (antiparallel to the vortex core) values. We 
demonstrate that in a certain range of the bias field mag- 
nitudes there are two stable gyrotropic modes of the vor- 
tex core rotation having different frequencies and oppo- 
site circular polarizations, and corresponding to opposite 
orientations of the vortex core relative to the direction of 
the bias magnetic field. The difference in frequencies of 



these two stable gyrotropic modes is proportional to the 
magnitude of the applied bias field, H, and, also, to the 
dot aspect ratio (3. We believe that this effect might be 
important for the development of novel magnetic memory 
elements. It allows one to determine the polarity of the 
vortex core by measuring the frequency of the resonance 
absorption in a nano-disk subjected to a perpendicular 
field, which is substantially easier to implement than the 
detection of the sense of signal circular polarization, as 
it was suggested in [1] . 

Our experiments were performed at room temperature 
on individual nano-disks of the thickness t = 43.8 nm 
and two different radii = 130 and i?2 = 520 nm). 
The nano-disks were patterned from a film of NiMnSb, a 
soft conductive magnetic material having very low mag- 
netic losses (typical Gilbert damping constant is a = 
0.002 - 0.003), cpitaxially grown on an InP(OOl) sub- 
strate A reference sample of a continuous film was 
cut out of the same film for characterization purposes. A 
50 nm thick SisN4 cap layer was deposited on top of the 
disks for protection, and a 300 nm thick Au broadband 
strip-line microwave antenna was evaporated on top of 
the disks. This antenna generates a linearly polarized mi- 
crowave magnetic field oriented in-plane perpendicularly 
to the stripe direction. The disk samples were placed in 
the uniform external bias magnetic field oriented perpen- 
dicular to the disk plane and having magnitude that was 
continuously varied from —10 to +10 kOc. 

The detection scheme of MRFM is inspired by mag- 
netic force microscopy (MFM) Q . It consists of an ultra- 
soft cantilever with a 800 nm diameter sphere of amor- 
phous Fe (with 3% Si) glued to its apex. The magnetiza- 
tion curve of this probe is typical of a soft Fe sphere (co- 
ercitivity < 10 G). Its magnetic moment always follows 
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TABLE I: Physical parameters of the reference NiMnSb layer. 
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FIG. 1: (Color online) (a) The cantilever frequency, propor- 
tional to M z , the vertical component of the static magneti- 
zation of the Tii = 130 nm NiMnSb nano-disk, as a function 
of the perpendicular bias magnetic field H ; (inset) Detection 
scheme: the vortex gyrotropic mode is sensed through the 
dipolar force induced on the spherical Fe probe; (b) Nano- 
disk excitation spectrum in the saturated regime (\H\ > H s ) 
recorded at 7.5 GHz; (c) Excitation spectra in the vortex (un- 
saturated) state (\H\ < H s ). 



the direction of the applied field if and senses a dipo- 
lar force proportional to the perpendicular component 
M z of the magnetization of the nano-disk. MRFM spec- 
troscopy is achieved by placing the center of the sphere 
above the center of the nano-disk. A ferromagnetic reso- 
nance (FMR) spectrum is obtained by recording the vi- 
bration amplitude of the cantilever as a function of if at 
constant microwave excitation frequency, that is switched 
ON and OFF at the cantilever resonance frequency. The 
MRFM signal originates from the diminution of M z of the 
nano-disk produced by the absorbtion of the microwave 
field Q. During the scan, the cantilever resonance fre- 
quency is also recorded (see FigQJi). This enables static 
magnetometry of the sample. The saturation field of the 
Ri = 130 nm nano-disk is H s = ±6 kOe, which marks 
the region of the maximum radial susceptibility. 

FigQJ> shows the FMR spectrum of the smaller nano- 
disk at the frequency of 7.5 GHz. Above the saturation 
field ff s , the excitation spectrum of the disk consists of 
a series of peaks corresponding to the confined dipolc- 
exchange spin-wave modes of the disk. Their quan- 
tized resonance frequencies (resonance fields) are given 
by Eq.(l) from jToj. These modes are labeled with a pair 
(I, m) of azimutal and radial mode indices . The lowest 
(0, 0) (and the most spatially uniform) spin wave mode 
of the disk is situated at 8.1 kOe and is marked by a blue 



circle symbol. The frequency of this lowest quasi-uniform 
Kittel mode ujk can be approximately described as: 



u K {H)= 1 {H-H K }, 



(1) 



where H K = 47rAi s (N zz - N rr ) - if a- Here M s is 
the saturation magnetization, Nu are the diagonal ele- 
ments of the effective demagnetization tensor of the uni- 
formly magnetized disk, Ha is a uniaxial perpendicular 
anisotropy field (easy-plane for Ha < 0), and 7 is the 
gyromagnetic ratio. All the parameters used in Eq. ^ 
were measured independently using the reference film 
sample by means of standard magnetometry and cavity- 
FMR techniques (results are presented in table fl} . The 
calculated position of the Kittel mode (as well as the 
positions of the higher modes calculated using Eq. (1) 
from [l(|) are shown near the corresponding index pairs 
in FigQJ). This calculation did not use any adjustable 
parameters and took into account the stray field of the 
cantilever probe (around 500 Oe) Q. The analytically 
calculated values of Hk for the smaller and larger disks 
are Hk ~ 5.7 kOe and Hk ~ 8.0 kOe, respectively. 

When the external magnetic field H is reduced be- 
low the saturation field H s , a single magnetic vortex is 
formed inside the nano-disk (which can be visualized by 
standard MFM). FigQJ: shows the dynamical behavior 
of the vortex-state disk (unsaturated regime). The ab- 
sorption peak observed in the spectra corresponds to the 
excitation of the gyrotropic mode in the disk (see inset 
of FigQJi) [ll[. It produces a mechanical signal mainly 
because the contribution of the vortex core to M z di- 
minishes when it is excited Q. Such signal is sensitive 
to the polarity of the vortex core: it is positive when 
the sphere magnetization and the vortex core are par- 
allel and negative in the opposite case. Moreover, due 
to the field dependence of the magnetic moment of the 
sphere, the amplitude of the MRFM signal decreases as 
I if I is reduced to zero. 

The baseline of each spectrum in FigQJ is set at the 
excitation frequency. The spectrum for 1.15 GHz has 
an intense peak at around +3.5 kOc (blue circle) corre- 
sponding to the excitation of the gyrotropic mode of the 
vortex core, as well as its mirror image at —3.5 kOc. The 
gyrotropic frequency decreases with the bias field if, and 
in the zero field (if = 0) has the value of 0.7 GHz. In the 
region of negative bias fields if < the mode continues 
to exist, but the corresponding resonance peak has a neg- 
ative amplitude, which indicates that now the direction 
of the bias field and the amorphous sphere magnetization 
is opposite to that of the vortex core. The gyrotropic fre- 
quency continues to decrease with the same slope until 
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FIG. 2: (Color online) Frequency of the lowest excitation 
mode in a magnetic disk as a function of the perpendicu- 
lar bias magnetic field for two sizes: (a) R\ = 130 nm, (b) 
7?2 = 520 nm. Dark blue (light red) dots are experimen- 
tal points corresponding to the field variation starting from 
large positive (negative) values. Solid lines are theoretical 
predictions from Eq. {1} for \H\ > H a and from Eq. ([2]) for 
\H\ < H s . 



the bias field value of H r = —2.8 kOe, where the res- 
onance peak abruptly changes its sign, frequency, and 
slope, which indicates reversal of the vortex core polar- 
ity. Below the switching field H < H r the gyrotropic 
frequency of the inverted vortex core increases like it did 
in positive bias fields before the core reversal. It is im- 
portant to note, that upon increasing the bias field, the 
gyrotropic mode persists at reduced frequency until the 
symmetric reversal field of H r = +2.8 kOe is reached, 
and the new reversal of the vortex core polarity takes 
place. 

The complete phase diagrams demonstrating the dy- 
namic hysteresis loops of the gyrotropic mode frequency 
in a perpendicular bias magnetic field for two disk sizes 
are presented in FigO The dark blue dots correspond to 
the case when the bias field is reduced starting from large 
positive values, while the light red dots correspond to the 
opposite case, when the field is varied starting from large 
negative values. 

Several important features emerge in these phase dia- 
grams: (i) For \H\ > H s the excitation frequency of the 
lowest mode in the disk is well described by the Kittel's 
Eq. {T]), and at \H\ = H s the field slope of the mode 
frequency changes abruptly, while the expected disconti- 
nuity of the mode frequency is not seen in the experiment; 
(ii) In the interval \H\ < \H r \ each value of the bias field 
(except H = 0) corresponds to two different gyrotropic 
mode frequencies corresponding to two opposite orien- 



tations of the vortex core relative to the bias field, and 
the difference between these frequencies increases with 
the bias field H and with the disk aspect ratio /?; (iii) 
The magnitude of the bias field H = H r , at which the 
reversal of the vortex core polarity occurs, seems to be 
approximately the same for both studied disks. 

The frequency of the gyrotropic mode in the external 
perpendicular bias magnetic field H can be calculated 
as a ratio ujg(H) = k{H)/G{H) of the field-dependent 
vortex stiffness n(H) to the field-dependent magnitude 
G(H) of the z-component of the vortex gyrovector using 
the method of Thiele's equation that was used in to 
calculate this frequency in the case when H = 0. Accord- 
ing to the general definition given in 01 (seeEq. (2) in 0) 
when the perpendicular bias field is introduced we can 
write the vortex gyrovector as G{H) = G(0)(1 —pcosO), 
where p = ±1 is the vortex core polarity, G(0) is the 
gyrovector in zero bias field defined in [7J], and 9 is the 
polar angle of the static magnetization at the disk lat- 
eral boundary. For a sufficiently large disk radius (much 
larger than the radius of the vortex core) is is possible 
to estimate 9 from the standard clcctrodynamic bound- 
ary conditions at the vertical boundaries of the disk as 
cos 9 = H/H s . 

The field-dependent vortex stiffness n(H ) is found from 
the assumption that the main contribution to the vortex 
energy comes from the dipolar interaction of the volume 
magnetostatic charges created by the in-plane magneti- 
zation component of the shifted vortex outside the core 
■ In this region the in-plane magnetization components 
depend of the direction of the external bias field as sin (9, 
and the stiffness n(H), that is proportional to the square 
of the radial derivative of the in-plane magnetization, is 
expressed as n (H) = n (0) sin 2 9, where k(0) is the vor- 
tex stiffness in the zero bias field defined in |7i]. Thus, we 
can write the explicit expression for the gyrotropic mode 
frequency in the perpendicular bias field as 



u> G (H)=uj G (0){l+p 



H 



(2) 



where u>g(0) is the gyrotropic mode frequency in the zero 
bias field defined by Eq. (7) of \v^. The approximate ex- 
pression for u>q(0) obtained in the limit of small aspect 
ratios [3 < 1 has the form uj g (0) = 5/(9tt)4ttiM s (3 0. 
The value of the saturation field H s can also be deter- 
mined as the crossing point of the dispersion curves (fTj) 
and J2J). It occurs at H s ~ Hk + 2ujg(0)/i found to be 
w 6.1 kOe for the smaller disk, a result that is close to 
the experimentally measured value shown in Fig[2^,. 

The theoretical curves calculated using Eq. @ are 
shown by solid lines in Fig|2l and it is clear from this 
figure that Eq. ([2]) gives a good quantitative description 
of the experimental data for two different values of the 
disk aspect ratio. This simple analytic theory does not 
account for the reversal of the vortex core polarity ob- 
served experimentally at H = H r (see Fig.©). There- 
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FIG. 3: (Color online) (a) Numerically calculated field depen- 
dence of the static magnetization of the NiMnSb disk (thick- 
ness t = 43.75 nm, radius R = 100 nm): the hysteretic be- 
havior of the static magnetization, M z , averaged over the disk 
volume (dashed loop) and over the vortex core (solid loop). 
The inserts show the spatial distribution of M z at four pro- 
gressively decreasing values of H. (b) Resonance locus cal- 
culated by the 3D micromagnetic code (dots) along with the 
analytical model (solid lines). 



fore, we performed micromagnetic simulations to repro- 
duce this experimentally observed effect. The results of 
these simulations are presented in Fig|3l 

Our micromagnetic numerical code [13[ calculates the 
stable configuration of the magnetization vector in the 
nano-disk by solving the Landau-Lifshitz (LL) equation 
in the time domain. A 3D mesh with a cubic cell of 
size 3.125 nm gives a discrete representation of the nano- 
disk magnetization. The nano-disk dimensions were cho- 
sen as R = 100 nm, t = 43.75 nm, while the values of 
the magnetic parameters were taken from table [I] The 
numerically calculated static hysteresis loop of M z (H) 
(where the vertical component of the static magnetiza- 
tion M z was averaged over the disk volume) is shown in 
FigJH^, by a dashed line. Color images of M z distribution 
across the radial section of the disk taken at four points 
corresponding to progressively decreasing values of the 
perpendicular bias magnetic field H are shown as inserts 
in Fig[3^i. Although the vortex core switching is not seen 
in this averaged static hysteresis loop, it is clearly re- 
vealed when the averaging of the static M z component is 
done over the spatial region occupied by the vortex core 
(~ 15 nm), as demonstrated by the solid lines in FigJ3^i, 
where the core polarity switching appears as a jump at 
\H r \ = 3.25 kOe. 

It should be noted, that the vortex core reversal in 
the negative perpendicular bias field was studied analyt- 



ically and numerically in 14J and was attributed to the 
formation of a Bloch point. Experimentally this effect 
was studied on an array of magnetic dots in [11. The 
magnitudes of the core reversal fields obtained in 14lll5| 
are consistent with corresponding magnitudes obtained 
in our experiment Figf2]and numerical modeling Fig[3^i. 

The dynamic hysteresis loop (i.e. the dependence of 
the gyrotropic mode frequency on the perpendicular bias 
magnetic field H) shown in Fig|3}3 was numerically cal- 
culated using a second code developed by the authors, 
which computes the full dynamic susceptibility tensor x" 
from linearization of the LL equation around the equi- 
librium configuration It clearly demonstrates the 
effects of the bias-field-induced dynamic bistability [l6| 
in the field interval \H\ < H r . Additional numerical cal- 
culations are in progress to study the details of the gy- 
rotropic frequency behavior in the vicinity of the points 
H = ±U S , where the vortex is formed, and H = ±H r , 
where the reversal of the vortex core polarity takes place. 
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